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ABSTRACT 

We present an in-depth, multi-wavelength study of the radio galaxy PKS J0334-3900, 
which resides at the centre of Abell 3135. The spectro-polarimetric radio observa- 
tions are combined with spectroscopic optical and X-ray data to illustrate the use of 
Head- Tail radio galaxies in revealing properties of the intracluster medium (ICM). 
Australia Telescope Compact Array (ATCA) observations at 1.4, 2.5, 4.6 & 8.6 GHz 
are presented with a detailed analysis of the morphology and spectral indices which 
give physical parameters to constrain the dynamical history of the galaxy. Using these 
constraints we produce a simulation of PKS J0334-3900. We find that this particu- 
lar Head- Tail morphology can be induced via a combination of orbital motion due 
to a binary companion and the relative motion through the ICM. New Chandra 
images of A3135 are presented from which we obtain a cluster electron density of 
n e,o — (1.06±0.11) x 10~ 3 cm~ 3 , a global temperature of 2.4 Ig'^g keV and a lower limit 
to the radio jet power of PKS J0334-3900 of 1.6 xlO 44 erg/s. A new redshift analysis 
of the cluster from available spectroscopic data demonstrates A3135 to be comprised 
of galaxies with 0.058 ^ z < 0.066 and gives a new mean cluster redshift of 0.06228 
± 0.00015. We also uncovered a background subgroup between 0.066 < z < 0.070. 
Additionally, ATCA Stokes Q and U data of Abell 3135 were used to obtain rotation 
measure values along the line of sight to PKS J0334-3900. Using our simulation we 
are able to infer the distance between the jets along the line of sight to be 154 ± 16 
kpc, which when combined with the difference in rotation measure between the jets 
provides a novel new way to estimate the average magnetic field within a cluster. A 
lower limit to the cluster magnetic field was calculated to be 0.09 ± 0.03 /iG. From 
these results, we have shown that different techniques can be combined from observa- 
tions of Head- Tail galaxies to infer information on the cluster environment, showing 
them to be an important class of objects in next generation all sky surveys. 

Key words: galaxies: active, galaxies: clusters: general, galaxies: clusters: individual: 
A3135, galaxies: general, radio continuum: galaxies 



1 INTRODUCTION 

Abell 3135 (A3135) is a galaxy cluster located at the north- 
ern edge of the Horologium-Reticulum supercluster (HRS). 
It is a relatively unstudied system, with little work done to 
date on the cluster dynamics. The most distinguishing fea- 
ture of Abell 3135 is the spectacular Head- Tail radio galaxy, 
PKS J0334-3900, found close to the cluster centre. Head-Tail 
(HT) galaxies are radio galaxies that have their jets bent in 
the same direction with the core elliptical galaxy and bent 
jets resembling a head and tail. The morphology and radio 
power of HT galaxies put them in the transition region be- 
tween the Fanaroff-Riley classes I and II ( |Fanaroff fc Riley 
1974). 



HT galaxies are found mostly within or on the periph- 
ery of galaxy clusters (though they are not only found in 
clusters, Xizhen et al. ( 1995 1). They have been found in re- 
gions of galaxy clusters that are dense and turbulent (Mao 
et al.|2009 1, at a rate of about 1-2 per cluster in the local 



(z < 0.1) Universe. In the more distant Universe, tailed ra- 
dio galaxies are believed to act as signposts for over densities 
in large-scale structure, with HT galaxies being associated 
with galaxy clusters up to z ~ 1 ( Blanton et al.|2003 \. There 
is evidence that this association continues to the limits of 
both cluster and tailed galaxy detection, of z ~ 2 ( |Dehgha"n"| 
et al.|2011 i. 



Since the time of their discovery 40 years ago (Miley 



et al. 1972 1 the morphology of HT galaxies has been be- 
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lieved to be the result of environmental interactions, com- 
monly due to the surrounding intra-cluster medium (ICM) 
I Burns 1998 ) . In particular, the bent morphology is the re- 
sult of ram pressure from the HT's relative motion with 
respect to the ICM. This ram pressure can arise via the in- 



fall into the cluster potential well (Gunn & Gott 19721 or 
from winds associated with "weather" in ICM which could 
arise either from mergers or mass accretion driven activity. 
We can therefore use the morphology of an HT to study the 
nearby cluster environment. 

In recent years, the use of tailed radio galaxies as en- 
vironmental probes has gained momentum as a method for 



cluster detection (Blanton et al. 2003 Smolcic 2007 Mao 



et al.||2010|>, examining the dynamics of individual clusters 



(jPfrommer fc Jones|20iTj ), measuring the density (|Freeland 



Enfil in|2003| |Guidetti et aL[2 008 ) 



20111 in the 



et al.||2008 | and velocity flows (Douglass et al. 
ICM, and for probing cluster magnetic fields (Clarke et al 
200Tj |Eilek fc Owen||2002| |Johnston-Hollitt||2003[ |Vogt fc 



However, to date, no study has combined the morphol- 
ogy and linear polarisation of an HT galaxy to improve es- 
timates of cluster magnetic fields. We present a novel new 
method, where we combine knowledge about the morphol- 
ogy and linear polarisation of PKS J0334-3900 to probe the 
cluster magnetic field. From the morphology of PKS J0334- 
3900, a simple simulation is used to study how the mor- 
phology was generated by the environment. From this sim- 
ulation, we can estimate the distance between the jets of 
PKS J0334-3900 along the line of sight. The distance and 
intervening ICM between the jets are used to explain the dif- 
ference in average rotation measure between the jets. The 
net rotation measure between the jets provides an estimate 
of the cluster magnetic field. 

The paper is presented as follows: Section 2 presents 
a spectroscopic analysis of A3135. Section 3 presents the 
current knowledge of PKS J0334-3900; Section 4 contains 
details of the radio and X-ray observations and the results 
found in this study of PKS J0334-3900; Sections 5 presents 
the simulation of PKS J0334-3900, made using the mor- 
phology obtained from the results, this simulation is then 
used with the rotation measure values to estimate the clus- 
ter magnetic field. Section 6 presents the conclusion. 

We adopt a standard set of cosmological parameters 
throughout with H = 73fcms~ 1 Mpc~ 1 , Q. m = 0.27 & 
Q A = 0.73. At the redshift of PKS J0334-3900, 0.062310 
± 0.000097 ( [Collins et al.|1995[ ), 1 arcsecond is 1.185 kpc. 
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Figure 1. Top: Histogram of redshifts for the 155 galaxies within 
one degree of A3135 and with redshifts between 0.058 and 0.075. 
The plot has been split into three subgroups with the following 
colours; white is 0.058 ^ z < 0.066, grey is 0.066 sg z < 0.070 
and black is 0.070 ^ z < 0.075. The dashed red line gives the 
redshift of PKS J0334-3900 showing it is in the central bin of 
the main cluster component (shown as the white subgroup). Bot- 
tom: scatter plot of available redshifts surrounding A3135, colour 
coded as per the binning in the histogram. The circle denotes the 
galaxies within one degree of the cluster which are included in 
the histogram above. Crosses show the location A3135 and other 
clusters in the HRS. 



2 OPTICAL PROPERTIES OF A3135 



Spectroscopic redshifts were collected from the literature to 
determine the cluster velocity distribution. From this distri- 
bution a new average redshift of the cluster was calculated, 
showing PKS J0334-3900 to be at the centre of this distri- 
bution. 



2.1 Spectroscopic Analysis of A3135 

There are 310 spectroscopic redshifts available in the liter- 
ature dGalli et al.||1993| iCollins et al 



2003||Way et al.|2005|pones|2009[ ) within 1 degree of A3135, 



1995| |Loveday||1996 
Shectman et al.| 19961 |Daiton et al.|1997 | |Katgert et al. 1998 
Ratcliffe et al.||1998||Tucker et al.[2000| Colless et al.||2001 



of these 155 have values between 0.058 and 0.075, which has 
been defined to be the extent of the "Kinematic Core" of 
the HRS ( jFleenor et al.||2005| ). Figure [l] shows the veloc- 
ity histogram of the 155 galaxies known to be associated 
with A3135 & the HRS (top panel) in addition to their re- 
spective position on the sky (bottom panel). The spectro- 
scopic data were split into three components, denoted on 
Figure[l]via different colours. Most of the galaxies have red- 
shifts between 0.058 and 0.066 (white bins) and mark the 
main part of A3135. The other two groups are located at 
0.066 ^ z < 0.070 (grey bins), which represent galaxies in 
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a nearby background group, and 0.070 ^ z < 0.075 (black 
bins) which are more distant galaxies separated by a signif- 
icant velocity gap, and thus are not expected to be part of 
A3135. As shown on the figure, a large number of the galax- 
ies with redshifts 0.058 ^ z < 0.066 are close to the centre 
of A3135, showing that the centre of A3135 is denser, as ex- 
pected. The over density of galaxies in the south-east is due 
to the nearby cluster Abell A3142, centred at RAj2ooo = 
03:36:17.3, DEC J2 ooo = -40:38:09 (marked by a cross on the 
plot). Also visible in the top of the lower panel of Figure [l] 
is A3145, denoted by a cross. All three clusters are members 
of the HRS. 

The most reliable published value of the redshift for 



A3135 is 0.0642 found using 82 galaxies (Struble & Rood 



19991. However, we find using 110 galaxies with redshifts 
between 0.058 ^ z < 0.066 (the white histogram in[l| that 
the average redshift is lower, with z = 0.06228 ± 0.00015. 
The higher previously published value of the redshift can be 
explained by the inclusion of the previously unknown sub- 
group at 0.066 ^ 2 < 0.070 in the Struble calculation. This 
new mean redshift is in much better agreement with the 
redshift of the central Dominant galaxy, PKS J0334-3900, 
shown as the red dashed line in Figure [I] falling within the 
peak cluster redshift bin. 



Table 1. Details of the observations of PKS J0334-3900 using 
the ATCA. Column 1 is the array configuration; column 2 is the 
date; column 3 is the integration time; column 4 gives the central 
observing frequency of the 128 MHz ATCA band and column 5 
is the resultant UV range listed as an interval. 



Conf. 


Date 


Time 


Freq. 


UV range 






(Min.) 


(MHz) 


(kA) 


6A 


2001 May 20 


59 


1384 


[0.5, 28.5] 








2496 


[1.8, 48.2] 






31.2 


4800 


[4,96] 








8640 


[6,172] 


1.5A 


2001 Aug 3 & 7 


76 


1384 


[0.3, 21.1] 








2496 


[0.6,36.2] 






63.1 


4800 


[0.8,23], [30,69.5] 








8640 


[3,43], [55, 128] 


750A 


2001 Jul 13 


79.7 


1384 


[0.1, 3. 5], [9, 17.9] 








2496 


[0.2,6], [15.4,30.5] 






55.1 


4800 


[0.5, 12], [30.2,59.5] 








8640 


[3,21] [55,105] 


375 


2001 Jun 15 


75.4 


1384 


[0, 2.1], [16.5, 28.5] 








2496 


[0.1, 4], [28, 49] 






64.5 


4800 


[0.4,8. 5], [56,96] 








8640 


[0.2, 13], [100,173] 



4 OBSERVATIONS AND RESULTS 



3 PKS J0334-3900 

Several different radio telescopes have been used to observe 
PKS J0334-3900, at frequencies ranging from 80 MHz to 8.6 
GHz. In particular, low resolution observations have been 
performed with the Very Large Array (VLA) at 1.4 GHz 



(Condon et al. 19981 & 4.89 GHz (Ekers et al. 19891 and 



the Molonglo Observatory Synthesis Telescope (MOST) at 
408 MHz JSchilizzi fc McAdam|1975] ) & 843 MHz ( |Jones fc| 
Mc Adam 111992 ). Previous radio observations have classified 



PKS J0334-3900 as a wide-ang le tail (WAT) galaxy , with a 
power consistent with an FRI ( jGovoni et al.|[2000[ ). In this 
work, we present higher resolution observations from the 
Australia Telescope Compact Array (ATCA) at 1.4, 2.5, 4.8 
& 8.6 GHz. 

PKS J0334-3900's optical host is a giant elliptical 
galaxy with an apparent Johnson blue magnitude of bj — 
16.00 ± 0.06 ( |Hambly et"aL]|2001a|b[). It has been identified 
as the brightest cluster galaxy ( Sun|200"9| , with a redshift of 
0.062310 ± 0.000097 ( [Collins et al.||1995[ ). The optical host 
has a redshift slightly less than the previous cluster average 



of 0.0642 ( |Struble fc Rood| [l999), but in good agreement 
with the new cluster redshift, 0.06228 ± 0.00015, derived 
here. 

PKS J0334-3900 is powered by a super massive black 
hole with a mass of 3.5 x 10 s solar masses (Bettoni et al 



20031, making it fairly typical for a radio galaxy (Woo & 



Urry 



2002). PKS J0334-3900 is at the centre of A3135, mean- 



ing it is in the most dense region of the cluster, this is con- 
sistent with HTs found in other clusters in the HRS (Mao 
et al.|2009 i. 



The results presented in this section were derived using the 
radio and X-ray observations of PKS J0334-3900 and the 
ICM respectively. The results derived from the radio & X- 
ray observations include: the Stokes I images, the spectral in- 
dex, polarisation angle maps, rotation measure values, and a 
Chandra X-ray image of the ICM. Each of these results con- 
tains useful information for understanding how PKS J0334- 
3900 interacts with its environment; the Stokes I images dis- 
play the morphology, the spectral index can be used to es- 
timate the age of the source, the polarisation angle maps 
help understand the morphology by showing the direction 
of flow within the jets, the rotation measures are required 
for understanding the cluster magnetic field, and the X-ray 
image is used to estimate the electron density. 

4.1 Radio Observations 

All ATCA observations of PKS J0334-3900 were extracted 
from the archive (PI - Saripalli). They were performed 
in 2001 using a phase centre of RA.12000 = 03:34:07.180, 
DECJ2000 = -39:00:03.190. The observations used the fre- 
quencies 1384, 2496, 4800 & 8640 MHz, and a bandwidth 
of 128 MHz over 32 channels. The channels on the edge of 
the band were removed, leaving an effective bandwidth of 
96 MHz. Due to cross-channel interference in the 1384 MHz 
observation, the 8 MHz interval centred at 1408 MHz was 
also removed. 

Each observation used four baseline configurations (see 
Table [T]), providing excellent uv-coverage. Taking all con- 
figurations into account, the total integration time was 290 
minutes at 1384 & 2496 MHz, and 214 minutes at 4800 & 
8640 MHz. 

The flux densities were scaled according to observations 
of PKS B1934-638, and the assumed total flux densities at 



© 2013 RAS, MNRAS 000,[T}{T7] 



4 Pratley et al. 




3 i 34 m 35 5 34 m 30 5 34 m 15 s 34 m 10 s 34 m 05 s 34 m 00 5 3 h 34 m 30 ! 34 m 10 5 34 m 00 ! 

Right Ascension (J2000) Right Ascension (J2000) 



Figure 2. Stokes-I 1384 MHz images of PKS J0334-3900. The image on the left is at full resolution (beam size 7.46" X 4.41") and has 
an RMS noise of 0.37mjy beam -1 . The image on the right has had a taper applied to give a beam size of 14.8" X 11.0" resulting in an 
RMS noise of 0.9 mjy beam -1 . In both images the contours start at 3 times the RMS noise and increase in multiples of \/2. 
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Figure 3. Stokes-I data at 2496 MHz for PKS J0334-3900 with a 
beam size of 4.1" X 2.5" and an RMS noise of 0.266 mjy beam -1 . 
The red contours show the shape of the jet from this Stokes-I 
image and increase in multiples of y/2 of the RMS noise until the 
6 th contour. 



1384, 2496, 4800 & 8640 MHz are 14.94±0.01, 11 14±0.01, 
5.829 ± 0.005 & 2.841 ± 0.006 Jy, respectively (|Reynolds 



19941. The time variations in complex antenna gains and 
bandpass were calibrated using observations of the unre- 
solved source PKS B0332-403. 

The observations were calibrated using MIRIAD, while 
following the standard ATCA pre-CABB reduction method 
(Sault et al. 1995[ ). The Stokes intensity images were also 



Table 2. Details of continuum Stokes I images made and shown 
in Figures[2][3]and^]and also used for spectral index calculations, 
but not shown. At 1384 and 8640 MHz two images were made, one 
at full resolution and one after applying a taper in the uv-planc. 



Frequency 


Synthesized Beam Size 


PA 


Stokes I a 


(MHz) 


(arcscconds) 


o 


(mjy b -1 ) 


1384 


7.459 x 4.406 


11.40 


0.367 




14.81 x 10.99 





0.867 


2496 


4.088 x 2.507 


9.826 


0.266 


4800 


1.799 x 1.434 


26.32 


0.150 


8640 


0.9669 x 0.7030 


9.132 


0.133 




2.617 x 1.996 


-79.15 


0.266 



restored using MIRIAD. A summary of the resolution and 



RMS noise for the Stokes I images shown in Figures [2j [3] and 
|4]can be found in Table [2] 



4.2 Radio Continuum Images 

The Stokes I continuum images show PKS J0334-3900 to be 
bent at an obtuse angle, the top jet is directed north and the 
bottom jet is directed south-east, each jet has a bright knot 
close to the core and ends with a clear asymmetric hook 
morphology (see Figure While there are similarities in 
shape between the jets, each jet varies differently in surface 
brightness. We find the jet's brightness variation is affected 
by its geometry, velocity, and orientation. For example, at 
1.4 GHz (see Figure|2|, the knot in the top jet is elongated in 
the jet's direction, but the knot in the bottom jet is brighter 
and is not as elongated. This can be explained if the bottom 
knot was directed towards the observer, it is expected to be 
compact in the plane of the sky and have a higher surface 
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Right Ascension (J3000) 

Figure 4. Digitized Sky Survey blue image of the core of PKS J0334-3900 and the bottom knot overlaid with 8640 MHz Stokes-I 
contours. Red contours are from a tapered image (beam size 2.617" X 1.996", RMS noise 0.266 mjy b -1 ), while blue contours give the 
highest resolution imaging and correspond to a beam size is 0.9669" X 0.7030" with RMS noise of 0.133 mjy b _1 . Contours start at 3 
times the RMS noise and increase in multiples of \f2. In each set of contours it is possible to observe a shape that looks like the bottom 
jet leaving the AGN. Note the presence of the close companion galaxy to the south west. 
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Figure 5. MOST 843 MHz Stokes-I map with a beam size of 
43" X 43" . Red contours arc the of the MOST 843 MHz Stokes I 
map, while the blue is the 3<r contour from the tapered 1384 MHz 
ATCA map presented in this work (see Figure [|| . At 843 MHz, 
there is more resolved low surface brightness emission at the ends 
of the jets, showing them to be longer than the 1384 MHz ATCA 
Stokes I image. 



brightness. This explanation is also consistent with the 2.5 
GHz observation (see Figure [3b, the knot in the top jet has 
the same elongation, however, the knot in the bottom jet 
has newly resolved structure. 

We now consider the surface brightness along the hooks 
and how this relates to their orientation. Figure [2] shows the 
hooks of the top and bottom jets. Each hook has a local 
peak in intensity, the bottom jet has a peak at the start 
of its hook, and the top jet has a peak half way around 
its hook. After the relative peaks, the bottom jet intensity 
decreases steadily around the bottom hook, and the top jet 



decreases to the end of the hook. For each jet, the intensity 
starts decreasing at different locations around their hooks. 
This difference could be explained if there is a difference in 
orientation between the hooks along the line of sight. If the 
end of the bottom hook is directed away from the observer, 
this could explain its steady decrease in intensity. This would 
mean that the bottom jet starts off directed towards the 
observer and changes direction when it reaches the bend 
around the hook. This explanation would suggest that the 
orientation of the top jet does not change as drastically with 
respect to the plane of the sky. 

Figure [4] shows the Digitized Sky Survey (DSS) image 
overlaid with Stokes I 8.6 GHz contours. At 8.6 GHz, the 
start of the bottom jet is detected at the radio core. The 
lack of a counter jet and the angle the bottom jet leaves 
the core suggests that the bottom jet is in the foreground 
(Parma et al. 1994). This detection is not visible at any 
of the other frequencies due to lack of resolution. Next to 
the core elliptical galaxy is a nearby companion galaxy (see 
FigureQ, with an apparent Johnson blue magnitude of bj = 
17.41 ±0.08 ( |Hambly et al.|2001a|b"| ). This companion has a 
projected distance of 20.54 kpc from PKS J0334-3900. 

Radio observatio ns of PKS J0334-3 900 by MOST (843 
MHz) (see Figure [S| ( |Bock et al.|[T999} and the NVSS (1.4 
GHz) show a similar morphology, while being at lower reso- 
lution. Unsurprisingly, the MOST observation has more dif- 
fuse low surface-brightness emission which cannot be seen 
in the 1384 MHz ATCA observation (see Figure pjj, this is 
likely due to the spectral steepening of the source. Since the 
MOST observation had more low surface brightness emis- 
sion, we used it to measure the projected linear size of each 
jet (see Table [3b. 

The physical attributes of PKS J0334-3900 were mea- 
sured and can be found in Table [3] this includes values of 
spectral index, flux, and total power. The uncertainties in 
the flux were calculated using AS = a 2 + (bS) 2 , where a 
is the noise of the image, b is the residual calibration error 



(assumed to be 0.01 ( Venturi et al.|20"00 l), and S is the flux 
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Right Ascension (J3000) 

Figure 6. Chandra observation of A3135 which has been adap- 
tively smoothed. The image has been corrected for background 
counts and exposure, and point sources removed. The blue con- 
tours are at the values 0.018, 0.024, 0.030, 0.036, 0.042, 0.048 & 
0.054 and the red shows the 1.4 GHz radio image of PKS J0334- 
3900. The red circle has a radius of 5.5', within which we have 
a reasonable kT temperature. The calculated estimate of n e was 
determined inside the purple circle of radius 7'. 



At 1384 MHz, PKS J0334-3900 has a total power of 
1.46 x 10 25 W Hz -1 , this places it at the transition phase 



between being an FR I and II source (Ledlow & Owen 1996 1, 
which is typical for HT galaxies. 



4.3 X-ray Analysis 

The X-ray observation of Abell 3135 was performed with 
the Advanced CCD Imaging Spectrometer (ACIS) on 2009 
May 7 (obsID: 9393). Standard Chandra data analysis was 
performed which includes the corrections for the slow gain 
change, charge transfer inefficiency and the ACIS optical 
blocking filter contamination. The observation was filtered 
for background flares and none were found. The total ex- 
posure is 15.3 ks. We used CIA04.4.1 for the data analy- 
sis. The calibration files used correspond to Chandra cali- 
bration database (CALDB) 4.5.0 from the Chandra X-ray 
Center. The solar photospheric abundance table by | AndeTs] 
& Grevesse ( 1989 1 is used in the spectral fits. Uncertainties 



quoted in the X-ray analysis are la. The 0.7 - 2 keV Chandra 
image of A3135, with background subtracted and exposure 
corrected, is shown in Figure [6] X-ray point sources have 
also been removed and filled with the average of the sur- 
roundings. The shown image was adaptively smoothed with 
the Science Analysis System (SAS) tool ASMOOTH. De- 
fault parameters of ASMOOTH Q were used, which includes 
a maximum width of the gaussian kernel of 19.68" (or 40 
ACIS pixels) and a desired signal-to-noise ratio of 10. 

For the spectral analysis, we use a Galactic absorption 



http: / / xmm.esac.esa.int / sas / current / doc / asmooth/index.html 



Figure 7. Three colour (blue, red, and infrared) image of A3135 
overlaid with Chandra observations (purple) and 1.4 GHz ATCA 
observations. The image colour scales have been adjusted in a 
similar manner to images made for the Hubble Heritage Project. 



column density of 1.6 x 10 2U cm~ J (Kalberla et al. 2005) 



The global spectrum is extracted within a radius of 5.5' 
and the point sources are excluded. We followed the method 
described in |Sun et al. ( 20091) for background subtraction, 



which involves the use of the ACIS stowed background and 
the modeling of the cosmic X-ray background. After the re- 
moval of the particle background and the cosmic X-ray back- 
ground, there are 1349 counts (0.5 - 7.0 keV) from the cluster 
in the global spectrum. The spectral fitting was performed 
with XSPEC 12.7.1 and the spectral model is APEC. The 
measured global temperature and abundance are 2.41 Ip'lg 



keV and 0.08 



solar respectively for AtomDB 1.3.1. The 



results with AtomDB 2.0.1 are very similar. C-statistic does 
not provide a goodness-of-fit measure but the XSPEC ver- 
sion of the C-statistic is defined to approach \ 2 m the case 
of large number of counts. The C-statistic of the fit is quoted 
for reference, 62.0 for 53 degrees of freedom. 

The radial surface brightness profile of the cluster, cen- 
tered on the central galaxy, is also measured in the 0.7 - 2.0 
keV band. The profile is analyzed with the CIAO/Sherpa 
tool and it is well fitted by a P model. Combined with the 
results of the global spectral fitting, the radial electron den- 
sity profile of the cluster is derived: 

n e = n e ,o[l + (r/r s ) 2 ]~ 1,5,s , 
n e , = (1.06 ±0.11) x 10~ 3 cm" 3 , 
13 = 0.504 ±0.016, 
r 3 = 133 ± 27kpc. 

Figure [6] shows the cluster X-ray emission is skewed 
to the southwest. There are several large elliptical galax- 
ies around the cluster center, including one (AM 0332-391 
NED02, RAjaooo = 03:34:06.7, DECj 2 ooo = -38:59:33) that 
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Table 3. Details of measured total length along the jets, flux and spectral index for the core, top jet, bottom jet and over the whole 
HT galaxy. The length was measured using the MOST 843 MHz observation. Linear least squares fits have been applied to determine 
spectral indices. Literature values used for the entire galaxy are shown in TableRland the core values are drawn from (|SIee et al.|1994|. 



Jet 


Linear Size 
(kpc) 


Frequencies 
MHz 


Flux 
(mjy) 


Polarized Flux 
(mJy) 


a This work 


a Inc. literature 


P 

Wffz" 1 


Total 


564 


1384 
2496 
4800 
8640 


1913 ± 22 
1368 ± 16 
937 ± 12 
445 ± 5 


280 ± 3 
190 ± 1 
28 ±0.2 


—0.79 ± 0.01 


—0.72 ± 0.01 


1.46 x 10 25 
1.05 x 10 25 
7.17 x 10 24 
3.40 x 10 24 


Top 


263 


1384 
2496 
4800 
8640 


785 ± 9 
697 ± 8 
438 ± 5 
222 ± 2 


90 ± 0.9 
100 ± 1 
10 ±0.1 


-0.71 ± 0.01 






Bottom 


301 


1384 
2496 
4800 
8640 


1120 ± 13 
841 ± 9 
488 ± 7 
210 ± 2 


190 ± 2 
90 ±0.1 
14 ±0.1 


-0.88 ±0.01 






Core 




1384 
2496 
4800 
8640 


8 ±0.8 

11 ± 0.3 

12 ±0.2 

13 ±0.2 




-0.15 ± 0.02 


-0.15 ± 0.02 





is ~0.13 mag brighter than PKS J0334-3900 at the K s band. 
It is likely that the system is in a late merging stage and the 
merger happened along the SW-NE direction on the plane 
of sky, which is consistent with the elongated optical distri- 
bution shown in Figure [l] 

Previous X-ray analysis has shown PKS J0334-3900 is 
associated with a small X-ray cool core (or a corona), (see 
Sun (20091, Table 2). There are no significant cavities de- 



tected in the X-ray emission. However, we can use the vol- 
ume of the radio lobes to estimate the cavity power, e.g. 
as in Birzan et al. (2008). If we assume a prolate spheroid 



for each lobe with 140 kpc x 40 kpc x 40 kpc for both the 
upper lobe and the lower lobe, using the average values for 
n e (6 x lO^cm" 3 ) and kT (2.4 keV), we find 4PV ~ 10 60 
erg. Assuming the sound speed for the bubble rising speed 
(McNamara & Nulsen 20071, the cavity power is then ~ 1.6 
x 10 44 erg s"\ which is reasonable for PKS J0334-3900's ra- 
dio luminosity. This value can also be taken as a lower limit 
for the jet power. 

We present the Chandra data in combination with both 
radio and a three colour (blue, red, and infrared) image in 
Figure [7] This image is presented to give a larger scale view 
of all emission in this system and provide an aesthetically 
pleasing example of the complexities of cluster emission. 



4.4 Spectral Indices 

Generally, the synchrotron spectrum of radio galaxies fol- 
lows the power law S oc i/ a , where S is flux, v is frequency, 
and a is the spectral index. The spectral index is useful as 
it determines the synchrotron spectrum for the source, and 
when used with the critical frequency, can be used to esti- 
mate the age of the synchrotron emission. 

Flux values for PKS J0334-3900 were measured from 
this work and collected from the literature, in order to deter- 
mine the spectral index. These values were plotted as \ogS 
vs log v, where the gradient is the spectral index. The total 



spectral index of the source from the ATCA data (see Figure 
[8| has a value of a — -0.79 ± 0.01. The literature has a large 
number of flux densities recorded for PKS J0334-3900 across 
an extremely wide range from 80 MHz to 8.6 GHz (see Table 
EJ. These literature values were combined with the measured 
values in this work (see Figure [8]), and the spectral index of 
the source was calculated to be —0.72 ± 0.01. 

Not all of the literature values were used in this calcu- 
lation, due to three factors. First, the flux density for 160 
MHz has been left out of this calculation, since the low reso- 
lution suggests nearby sources are blended with PKS J0334- 
3900, making it unreliable as a flux estimate. Second, where 
there are multiple flux densities at the same frequency, we 
selected the most reliable value for the calculation (see Ta- 
ble EJ). Finally, the 178 MHz data point was not included 
in the calculation (seen as red in Figure |8|, because it is an 
interpolated value. 

Figure [8] shows that the point at 80 MHz appears to 
show some deviation from the power law suggesting the crit- 
ical frequency to be near the interpolated point at 178 MHz. 
If we assume 178 MHz is the critical frequency, we may use 



the method deployed in Lai & Rao (20041 to calculate an 



estimated upper limit age of the synchrotron emission using 
equation (1) (the details behind this calculation are given in 



Miley (19801) 



t = 1060- 



B i/a 



(1 + 2) 



-1/2 



I'br 



-1/2 



(1) 



where t is the upper limit to the age of the jets in Myr. 
B and Bic are the total and inverse Compton derived mag- 
netic fields respectively, both measured in fiG, and V\, r is 
the break frequency of the injected electron population in 
the jets. 

This method assumes the magnetic field magnitude for 
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log(S) = -0.721og(y) +12.72 
Aog(S) = -0.791og(y) +13.32 



lu R (,/> 



Figure 8. Total spectral index plot of PKS J0334-3900 using only 
the ATCA data from this work (green), and with all available 
values from literature (blue). The equation of the solid trend-line 
for the ATCA only data is log 5 = -0.79 logy + 13.32, while 
the ATCA and literature values produce the dashed trend-line, 
with the equation logS 1 = -0.72 logy + 12.72. The ATCA only 
line corresponds to data in Table [3] whereas the literature values 
along with the ATCA values from this work can be found in Table 
[4] Note that the red data point represents an interpolated value 
from literature and was not used in the line fitting. 



B = 



(2) 



the inverse Compton component is given by 

Bic 

and is a function of redshift and can be estimated as Bi C — 
4(1 + zf for the case of maximum radiative lifetime of the 
particles. These equations yield a magnetic field value of 
B = 2.6/iG and an upper limit age of 146 Myr for the jet 
emission. 



4.5 Polarisation 

Stokes Q and U images of PKS J0334-3900 were made from 
the 1384 and 2496 MHz observations, and used to generate 
maps of the linear polarisation angle (PA) and the frac- 
tional polarisation (see Figures [9| |10] & and Table ^. 
These maps are useful for understanding the jet morphol- 
ogy. The PA maps have their angles rotated by 90° to show 
the direction of magnetic field along the jets. 

Figures[9]&[l0]show the magnetic field (B) vectors along 
the bottom and top jets respectively. The PA values are 
masked if their error was greater than 10° , or if the polarised 
intensity was below 3 sigma. At each frequency the vectors 
point along the jets. 

At the start of the bottom jet (see Figure Rjb, the 2.5 
GHz observations show the magnetic field leaving the knot, 
and curving steeply around the start of the hook. This steep 
curvature of the magnetic field corresponds to where there is 
a peak in intensity. This may suggest that the jet is directed 
towards the observer just before curving into the hook. Just 
after the hook starts, there is a break in the polarisation 
signal at RA.J2000 = 03:34:14. The break may be a sign of 
the jet changing direction. At 1.4 GHz, the vectors are per- 
pendicular to the break at the start of the hook, and then 



Table 4. All frequencies and net flux densities of PKS J0334-3900 
currently known, along with ones measured in this paper. * this is 
not a measured value, but rather is an interpolated value between 
literature low frequency measurements. We list the interpolated 
value here for completeness, but do not use it in calculations of 
the total source spectral index. We have assumed a 10 per cent 
error for this interpolated source, o claimed in the literature to be 
from |Schilizzi fe"M cAdam ( 1975), the correct value is not found in 
that paper but rather in |Wa ll & Schilizzi (1979). * measured from 
the Sydney University Molonglo Sky Survey image, t measured 
from the NVSS image. There is a corresponding graph, Fig. [8] 
that includes all the frequencies denoted with a 1 , but does not 
include the flux density at 160 MHz as this was believed to be 
unreliable. 

Frequency Flux Density Reference 



MHz 


mjy 


80 1 


10000 ± 1500 


160 


3700 ± 555 


178 


7300 


408 1 


4200 ± 210 


635 1 


3820 ± 96 


843 


2410 


843 1 


2578 ± 25 


960 1 


2430 ± 48 


1384 1 


1913 ± 23 


1410 


1700 


1410 1 


1790 ± 60 


1420 1 


1570 ± 78 


1477 1 


1642 ± 16 


2496 1 


1368 ± 16 


2700 


900 


2700 1 


940 ± 28 


4800 1 


937 ± 12 


5000 1 


740 ± 22 


8640 1 


445 + 5 



Mills et al. (1960) 
Kkers et af] dl989 1 
[Burge ss & Hunstcad {2006 I 
I Wall fc ijcMlizzi| ( |1979p 
Wills||l975} 
I J ones &; McAd am 
This work* 



(1992) 



Wills 



1 1975 1 



This work 
|Bolton et al 
|wills|(l975 
|Fkers et al. 
This work' 
This work 



(1965) 



(1989) 



Bolton & Shimi 


nins 


(1973 


Wall &; Cannon 


{197 





Wall & Schilizzi 
This work 



1 1979 1 



follow the curve of the hook thereafter. However, at 2.5 
GHz, the vectors appear to be curving into the hook just 
before the break. If the break in polarisation is where the 
hook changes direction along the line of sight, it is consis- 
tent with the surface brightness decreasing around the hook, 
as explained in Section |4.2| At the end of the hook, there 
is a large, ordered field structure and high fractional po- 



larisation (see Figure 111. Additionally, the bottom jet has 
significantly higher fractional polarisation when compared 
to the top jet. The 2.5 GHz observation suggests that the 
end of the bottom jet could be curving down, however the 
1.4 GHz image shows that the situation is more complicated, 
suggesting upward motion. However, it is well known from 



the Laing-Garrington effect (Laing 1988 Garrington et al. 
|1988[ ) that the foreground jet in a radio galaxy will exhibit 
higher fractional polarisation, suggesting the bottom jet and 
in particular the end of the jet is directed towards the ob- 
servers, and thus is likely to be moving up and towards the 
line of sight. 

The top jet has a curved flow out of the knot and bends 



around the hook of the top jet (see Figure 10 1. The intensity 
image does not show the jet curving out of the knot, how- 
ever, it is consistent with the observed hook shape. At the 
end of the hook shape, the 1.4 GHz vectors suggest that the 
jet flows down at a 45 degree angle. The lack of PA vectors 
at the end of the jet make it difficult to interpret the flow. 
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Figure 9. The direction of the magnetic component of the polarised emission is compared to the morphology of the bottom jet. The 
lefthand images are at 1384 MHz, and the righthand images are at 2496 MHz. The contours on the intensity images start at 3<r then 
increase as multiples of y/2. The image details are given in Table [5] Note the downward 'finger' of emission in the 1384 MHz image at 
RAj2ooo 03:34:11.06 is due to the AGN associated with cluster member AM 0332-391 NED04, which is clearly seen as a discrete source 
in the 2496 MHz image on the righthand side. 




Figure 10. The direction of the magnetic component of the polarised emission is compared to the morphology of the top jet. The first 
two images on the left are at 1384 MHz, and the last two on the right are at 2496 MHz. The contours on the intensity images start at 
3<t then increase as multiples of \f2. The image details are given in Table pi 



4.6 Rotation Measures 

Faraday rotation is the phenomenon of linearly polarised 
light travelling along a magnetic field causing a rotation in 
polarisation angle (%). This rotation depends on the electron 
density (n e ), magnetic field magnitude along the line of sight 
(B||), distance passed through the magnetic field, and the 
wavelength of the linearly polarised light (A). Rotation mea- 
sure, RM = jA, follows the relationship RM oc J n G B|| • df. 



The value of RM can be obtained by measuring \ for differ- 
ent values of A. As tailed radio galaxies present extended lin- 
early polarised structures, they have been used in studies of 
cluster magnetic fields using Faraday rotation ( Clarke et al 
[2001] [Euek fc Owen||2002| |Vogt fc Enfflin||2003| |Johnston- 
Guidetti et al.|2008[ ). ~ 



Hollitt 2003 



Using two different methods, the rotation measure 
across PKS J0334-3900 was calculated using the 1384 MHz 
Stokes Q and U observations. The first method was rotation 



© 2013 RAS, MNRAS 000,[T}{T7] 



10 Pratley et al. 

Table 5. Details of continuum Stokes I, Q and U images shown in Figures |lO| fc [9] 



Frequency 
(MHz) 


Synthesized Beam Size 
(arcseconds) 


PA 

o 


a(I) 

(mjybeam ) 


c(Q) 

(mjybeam ) 


«r(U) 

(mjybeam ) 


1384 

Polarised Image 
Intensity Image 


8.749 x 5.044 
15 x 15 


9.738 



1 


0.174 


0.122 


2496 

Polarised Image 
Intensity Image 


4.977 x 2.871 
15 x 15 


7.897 



0.833 


0.197 


0.192 




measure synthesis (Brentjens & de Bruyn 20051, and the 



second the traditional measurement of the gradient of the 
PA with frequency ( Burn|1966 1 . The RM values from these 
methods were then compared. Finally, a RM synthesis map 
of PKS J0334-3900 was made, and using the vector averaged 
polarisation an RM for each jet was calculated. 

The RM values were measured within a bandwidth of 
96 MHz, using 8 MHz intervals. Each interval had a contin- 
uum Stokes Q and U map which were combined into cubes. 
Only 12 intervals were used, because the interval centred 
at 1408 MHz was removed due to cross-channel harmonics. 
The Stokes Q and U maps (planes in the cube) have values 
of RMS noise between 330 & 380 njy, and the Stokes Q and 
U RMS values were added in quadrature to obtain the RMS 
noise in linear polarisaton. 

Values of RMuncar were calculated using the method 
of a least squares fit of \ vs - 

A 2 . A 2 was determined by the 
mean frequency of each interval. Values of RMuncar were 
masked when 5 or more of the points in A 2 -space had linear 
polarisations less than 3er. x values were corrected by nn 
where there was an obvious ambiguity. As the intervals are 
closely spaced, there were no instances where more than a 
single 7t correction needed to be applied. 

Values of RMsynthesis were measured using the RM syn- 
thesis technique, from the same 8 MHz intervals. These in- 
tervals were used to calculate the rotation measure transfer 



function (RMTF), shown in Figure [12] This RMTF has a 
full width at half maximum (FWHM) of 542.4 rad m -2 . The 
FWHM of the RMTF is large, so no attempt was made to 



deconvolve between components in RMsynthesis, because it is 
unlikely that the components would be resolved. Since mul- 
tiple components could not be resolved, any possible mag- 
netic fields in the foreground of the cluster cannot be distin- 
guished. The RMsynthesis values were measured by finding 
the location of the peak in the Faraday dispersion function. 
This was done by fitting a 2nd order polynomial to the peak 
over an interval of ±100 rad m -2 . These RM values were 
then used to generate an RMsynthesis map of PKS J0334- 
3900, shown in Figure ^3] 

To compare the two RM methods for each pixel value 
a scatter plot of RMLmear vs. RMsynthesis was made (see 

values are expected to yield su- 
-space with less 
This scatter plot 



Figure 141. The RMsynthe 



perior results to a straight linear fit in A 

[2012|. 



spurious results (Macquart et al 



conforms closely to the expected RM Linear — RMsynthesis re- 
lationship for a single source of RM along the observer's line 
of sight. This shows that the two methods are consistent on 
a pixel by pixel basis, as is expected for 97 per cent of sources 
with signal-to-noise ratio in excess of 3<r in Stokes Q and U 
(Macquart et al. 2012 I. Additionally, the close agreement be- 



tween the two methods demonstrates that the nir ambiguity 
correction applied to the RMLincar values was reasonable. 

The pixel by pixel RMs show good agreement with a 
A 2 fit in the linear plots (not shown) and have only a sin- 
gle, albeit wide, peak in the RM synthesis results. Addi- 
tionally, there is little depolarisation across the jets. Both 
the good A 2 fit and the lack of depolarisation suggest that 
the RM measured is not from the jets of the tailed radio 
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Figure 13. Masked rotation measure image of PKS J0334-3900 made using rotation measure synthesis. The mask corresponds to the 
same mask as on the polarization angle images in Figure [9] and Figure |To"| . The beam size of this rotation measure image is 7.1" X 7.1". 



ilar to that shown in Eilek & Owen ( 2002 1 , which presents a 



RM IlnI „ = 1.015(RM SvIltllalI )-1.162 



-150 -100 



-50 50 

RM Hyatlieaji (rad rrr 2 ) 



Figure 14. Comparison of measurements of rotation measure 
generated via rotation measure synthesis and linear fits per pixel 
for PKS J0334-3900. The good agreement between the two meth- 
ods is suggestive of a single component in RM space, and provides 
confidence that the methods were applied correctly. 



galaxy itself, but rather the tailed radio galaxy is providing 
polarised light which passes through a foreground Faraday 
screen. The reasoning for this is as follows: 1) if the RM was 
due to the jets, there would be different Faraday screens due 
to electron density and magnetic field variation through the 
interior of the jet, meaning there would not be a linear \ vs - 
A 2 fit. 2) if the RM was from the jets, we expect depolarisa- 
tion, as monochromatic waves from different depths into the 
jet would under go different amounts of Faraday rotation, 
meaning the angles would not be coherent. This result is sim- 



RM map of binary HT galaxies in Abell 400, in which a jet 
from each HT galaxy crosses over the other. From the RM 
maps, it is possible to see structure preserved between the 
jets seamlessly, suggesting the RM structure is not from jets, 
which should have independent RM structure, but rather 
from a foreground Faraday screen. Such a screen could ei- 
ther be from either the Milky Way or the cluster in which 
the sources are embedded. PKS J0334-3900 is far above the 
Galactic plane, with a latitude of —54.3°, so the magnetic 
field of the Milky Way is not expected to make a significant 
contribution to the RM. Estimations of the Galactic RM in 
the foreground of t he cluster give values o f 3.25 and 3.62 
± 5.65 rad m 2 , from 



Johnston-Hollitt et al. 



(20041 and Op 



permann et al. (20121 respectively, leading us to conclude 



that the observed RM is due to the magneto-ionic plasma 
in Abell 3135. We find the RM structure along each jet is 
coherent at a scale of order 15" which is two times the beam 
size of 7.1" x 7.1", as shown in Figure 13 The average and 
standard deviation of RMs yn thcsis for the jets can be found 
in Table H 

An RM was calculated from the vector average linear 
polarisation of the jets. To do this, the jets were unresolved 
by using a beam size of 250 arcseconds. Then the Stokes 
Q and U of each jet was measured to calculate the vector 
averaged P for each jet. From the two jets, the rotation 
measure was calculated using RM synthesis, using the same 
12 intervals. The resultant RM was found to be 27.4 & 5.5 
rad m -2 from the bottom and top jets respectively. Table 
[6] compares these RMs, and the average RMs from the RM 
synthesis map. There is a difference between the P-vector 
average RM and the simple pixel based statistical analysis. 
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Table 6. Table comparing the RM values found in the top and 
bottom jets, from statistics applied to the pixels, and the vector 
average polarisation, fi and a are the mean and standard variation 
for the rotation measure derived using the pixel based analysis, 
assuming Gaussian statistics. 



Jet 




a 


P-vector Average RM 




rad m — 2 


rad m -2 


rad m — 2 


Top 


1.16 


29.05 


5.5 


Bottom 


7.44 


24.88 


27.4 


Both 


5.54 


26.37 





In this case it was found that the P-vector average RM was 
dominated by a small number of high polarisation pixels. It 
is therefore more correct here to use the pixel based values to 
determine the average RM over each jet e.g. |Guidetti et al.| 
( 20081) &|Guidetti et all (12010 1. 



5 DISCUSSION 

5.1 Morphology of the Head- Tail galaxy 

PKS J0334-3900 has a morphology consisting of two bent 
jets and asymmetric hooks, which can be seen in Figure [2] 
Since both jets have a similar hook shape, the hooks were 
likely created by a mechanism which affects each jet simul- 
taneously. If the creation of the hooks was solely due to 
the motion of the ICM, the ICM would affect the jets in- 
dependently and asymmetrically. A simpler explanation for 
how the hooks were generated can be found by consider- 
ing a possible orbital motion of the AGN and its compan- 
ion (see Figure [4| . This case would intuitively give a solu- 
tion, which generates the similar scale and asymmetry be- 
tween the hooks. The presence of close companions to HTs 
has been confirmed spectroscopically in the literature ( |Rose| 



1982 Mao et al. 20091, and it has been suggested that a 
close companion is a necessary, but not sufficient condition 
to generate their morphology. If the possible orbital system 
here was falling into the cluster, or there were large-scale 
flows in the ICM due to a cluster merger, the HT galaxy 
would be affected as if there was a 'wind' in the ICM. This 
wind would cause both jets to bend as a whole in the same 
direction, preserving the hook shapes generated from orbital 
motion. 



5.2 Modelling of the tailed radio galaxy 

A simple simulation was performed to test if the hook and 
bent-tail morphology can be induced via orbital motion and 
a cluster wind. This model consists of a spinning and rotat- 
ing source emitting a steady stream of spheres into two oppo- 
site directions. We establish a coordinate system, x\, X2,Xs 
in which the model is centred on the AGN-companion mass 
centre, and defined such that the orbital plane lies on the 
ooi — X2 plane, and ipo — 1.67 rad sets the initial azimuthal 
position of the galaxy relative to X2 axis (see lower right side 
of Figure 15 I. Our model consists of symmetric and, to a 



Colina 1993 C larke fc Harris|1996 1. On kpc scales, the sym- 
metric structure and luminosity of the jets in FRI sources 
suggest that the jets quickly decelerate to non-relativistic 
speeds (Laing & Bridle 2002), this is supported by both the 
radio emission on large-scales and the X-ray analysis. 

The equations of motion of the spheres with the pro- 
duction rate of R can be written as 

Xlk = («{ fc + v° k + vf k )(t - ~) + x° lk 

for k — 1,2, and 3, where v 3 lk , v° k , and v]" k are initial velocity 
components of the I sphere due to jets, orbital motion, and 
wind respectively. These velocity components are defined by 
the following equations: 

v J cos(d) cos(4>) k = 1 
v 3 cos(#) sin(</>) k = 2 , 
v 3 sin(6») fc = 3 

-tp'r sin(i/)) k = 1 
v°k = { tp'r cos(^) k = 2 
k = 3 



and Vik = Constant, 



where 8, <f>, and ip are the polar and azimuthal angles of the 
precession axis, and the azimuthal angle of the AGN in the 
orbit respectively. x° k represents the initial position of the 
AGN in the orbit. 

In this simulation we assumed the jet's velocity to be 
~1 per cent of the speed of the light yielding 7 ~1, so the 
relativistic correction is not important. The model does not 
include energy loss or particle interactions. In Blandford & 
Icke ( 1978 1, a similar model was applied to the radio galaxy 



first approximation, ballistic jets moving at constant speed. 
Ballistic jets have been frequently used as a first approxi- 
mation for radio galaxies ( Blandford fc Icke||1978 Borne & 



3C 31, to test if the hook shapes in 3C 31 were due to or- 
bital motion with a companion. However, the model of 3C 
31 does not consider an interaction with the ICM, such as 
a 'wind', or consider the possibility of the radio galaxy un- 
dergoing precession. The model of PKS J0334-3900 includes 
such variables required to generate the hook and bent tail 
morphology. 

Parameters for the model were chosen using the avail- 
able radio, optical and X-ray data. The radio structure and 
polarization vectors of the extended radio source were used 
to match the morphology, and the visible luminosity of the 
sources were used to calculate their mass. The model's initial 
conditions and parameters were determined by ensuring the 
best possible fit to the mean width of the radio jets. These 
parameters are given in Table [7] and the model is displayed 
in Figure 15 (top). 

The major wide angle of PKS J0334-3900 emerges from 
its expanding tails' interacting with the surrounding ICM, 
or alternatively, it can be related to large-scale flows in the 
ICM. The relative motion of the HT through the intra- 
cluster medium generates strong ram pressure on the tails, 
giving the required angle between the jets. 

The difference in scale between the top and bottom jet 
cannot be explained by viewing angle alone. The ratio of 
the hook length to that of the total jet length is the same 
in each, despite the jets having different total lengths. This 
means the foreshortening of the top jet seen in the image 
cannot be produced by rotation of the viewing angle in the 
simulation solely. Furthermore, this can be explained if there 
is a difference in the actual length of the jets. This can be 
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Table 7. Orbital and characteristic parameters for the PKS 
J0334-3900 model. 



Observed Parameter 


Value 




Red magnitude of PKS J0334-3900 


14.92 ± 0.06 




Red magnitude of the companion 


15.92 ± 0.08 




Projected separation 


20.54 kpc 




Dynamical Parameter 


Value 




J0334-3900 mass 


5.23+2-qo X 10 12 


M Q 


The companion mass 


1.44±°;g x 10 12 


Mq 


Actual separation 


35.21 kpc 





Orbital velocity 194.88 km/s 

x\ component of radial velocity 516.54 km/s 

X2 component of radial velocity 1162.21 km/s 

x% component of radial velocity -86.09 km/s 

Jet velocity 2060.98 km/s 

Counterjet velocity 1613.53 km/s 

Eccentricity 

Initial orbital phase 95.73° 

obliquity 8.10° 

Precession period 84.0 Myr 

The dihedral angle 65.63° 

Evolution time 145.9 Myr 



described by variance in the jet and counter-jet velocity, or 
the ICM's density gradient; the top jet encounters a denser 
region. The latter is consistent with the overlay of the X- 
ray image shown in Section |4.3[ which shows the end of 
the foreground (bottom) jet lies outside the X-ray contours, 
which is not the case for the top jet. 

The minor features of the tail in the simulation are gen- 
erated by the close encounter of two galaxies; the host galaxy 
and the smaller companion galaxy (see Figure [4| orbit a 
common centre of mass. Mass estimation of the system was 
made following the assumption that there is a power law cor- 
respondence between the halo mass and the red luminosity 
of the galaxy: 



M = K'h 



with Mi, 



10.73 ± 2.53(10 11 /i" 1 M G 



1.51 ± 



0.04(10 10 /i~ 2 L©), and f3 = 1.4 ±0.2, optimized for the early 



type galaxies (Guzik & Seljak 20021. Unfortunately the r 



band rest frame luminosity of the galaxies could not be cal- 
culated via fitting of a template spectral energy distribution 
(SED) due to the lack of the other luminosity bands. How- 
ever, examination of k-corrections for other r-band magni- 
tudes at a redshift of 0.06 demonstrates at worst a maximum 
correction of 0.1 would be likely ( Blanton fc Roweis||2007 1, 
while a k-correction based on an empirical polynomial fit 
( jChilingarian et al.|20"To| |Westra et al.|2010| ) would result in 
something more like 0.08. As the uncertainty in the mass es- 
timate is dominated by other factors, the application of even 
the largest likely k-correction at this redshift gives a mass 
which is almost indistinguishable from a mass derived from 
an un-corrected r-band magnitude. Without a proper SED 
fit, these other methods introduce other assumptions, which 
may not hold. For this reasons, non rest frame red magni- 
tudes were used for the mass calculations which is given in 
Table [7] These estimations give rise to additional free pa- 
rameters, such as the dark matter portion of the general 



mass in the galaxies and also dynamical attributes such as 
eccentricity of the orbit. Nevertheless, uncertainty in these 
estimations would not create major issues since this is a 
qualitative way to understand the responsible mechanisms 
forming the morphology. 

It is essential to emphasize that the proposed model 
is a qualitative way to explain the morphology of this HT. 
In other words, due to a number of free parameters, which 
were not possible to determine from the available radio and 
optical data, this model is merely one possible way to justify 
the morphology. 



HT galaxies as probes of the intra-cluster 
magnetic field 



5.3 



Using the model in Section [5.2| we find the angle between 
the jets along the line of sight is 120 ±5° , the linear length is 
346 ± 32 kpc in the bottom jet and 170 ± 16 kpc in the top 
jet. So we find the distance between the jets, at halfway 
down them, is 154 ± 16 kpc. Using an average valu e of 
electron density, 1.09 x 



4.3 



10 cm derived in Section 
and the difference in RM between the jets, we find an av- 
erage intra-cluster magnetic field, along the line of sight, to 
be 0.05 ±0 02^G. Supposing that (B) = V3{B l} ) ( |Eilek fc 
|Owen|[2002[ ), we find the average magnetic field to be 0.09 
±0.03/xG. However, this is for the simplest case of a uniform 
magnetic field within the cluster, for more complicated fields 
it is expected to be larger due to random rotations in the 
cluster. Eilek & Owen (20021 shows how it would be possi- 



ble to compare the magnitude of B for different models of 
magnetic fields causing the Faraday rotation, however, that 
is beyond the scope of this paper. 



6 CONCLUSION 

In this work we present new radio, optical and X-ray anal- 
ysis for A3135 and the associated HT galaxy PKS J0334- 
3900. We present a new cluster mean redshift of 0.06228 ± 
0.00015, which is lower than the previously published value 
due to the presence of a heretofore undetected subgroup, 
with 0.066 ^ z < 0.070. We present new Chandra images of 
the clusters which reveal an extended X-ray morphology and 
provide evidence of a recent merger. Using new radio analy- 
sis at multiple frequencies we present a detailed analysis of 
the central HT galaxy PKS J0334-3900. We show that the 
morphology and use of linear polarisation of an HT galaxy 
can be used to determine not only the properties of the HT 
itself, but also as a probe of conditions in the ICM. The mor- 
phology and estimated age of the jets were used to make a 
simulation. This simulation showed that if the HT is falling 
into the cluster along with its companion about which it is 
undergoing an orbital motion, the observed morphology can 
be produced. Alternatively, large-scale flows in the ICM due 
to a cluster-cluster merger in combination with the orbital 
motion of PKS J0334-3900 and it's companion about their 
centre of mass can produce the same morphology. This sim- 
ulation was used to derive a likely angle between the jets 
along the line of sight to be 125 degrees. Furthermore, we 
estimate the upper limit on the age of the jets as 146 Myr, 
which also provides a time frame where the direction of the 
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145 . 9Myr 





Distance to the galaxy centre (kpc) 



50 100 150 200 250 300 350 400 450 




Figure 15. Top figure: Multi-view orthographic projections of the PKS-J0334-3900 model. The transparent region corresponds to 
the gap in the radio image close to the AGN. The top/bottom jets are moving away/towards us respectively. Lower left figure: 3D 
orientations of the orbital plane (O) and plane of sky (S). Arrows in the top right of the box show the direction of the presumptive wind. 
Lower right figure: Relative orientation of the coordinate system used in the model (x 1X2X3), and the equatorial coordinate system 
(xyz). The orbital plane lies on the x\ — X2 plane, mc and g represent the mass centre of the system and the radio galaxy respectively. 
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'wind' resultant from infall is recorded in the morphology of 
the jets. 

Additionally, rotation measure maps for each jet showed 
that the magnetic field was ordered on scales of roughly 15" 
(18 kpc), and the averaged RMs for each jet were 1.16 and 
7.44 rad m -2 for the top and bottom respectively. We at- 
tribute this difference to the different path lengths of light 
from each jet passing through the ICM. Using the linear 
lengths of the jets and the angle derived from the simula- 
tion, we obtained the likely distance between the jets along 
the line of sight. This distance was used in conjunction with 
the difference in average RM between the jets and average 
electron density, to calculate an average cluster magnetic 
field along the line of sight. This gave a lower limit of 0.05 
±0.02/iG for the line of sight component which suggests a 
lower limit for the total cluster field of 0.09 ±0.03^G. While 
the magnetic field is slightly weaker than typical literature 
values for cluster fields using other methods, which show 
ranges between 0.1 to a few fiG dependent on different loca- 
tions in the cluster (see for example comparison of different 
values across A3667 ( Johnston-Hollitt 2004 1), as our value 
is a lower limit, this does not present a problem. This is a 
novel method for estimating the cluster magnetic field which 
can be added to the repertoire of methods to probe cluster 
magnetic fields. 

In the next decade a number of radio telescopes 
are expected to be built or upgraded including; LOFAR 
( Rottgering et al.|2003| ), MWA ( |Tingay et al.|2013| ), ASKAP 
( Johnston et al.|2007[ ) ' MeerKAT~( |Jonas|2009[ >, WRST and 
the EVLA. These instruments are anticipated to conduct 
a series of all-sky continuum radio surveys such as EMU 



(Norris et al. 20111 and WODAN (Rottgering et al. 2011) 



and undertake a wide range of survey science projects (Nor- 
|ris et ah||2013[ ). These next generation surveys will produce 
catalogues of radio galaxies numbering in the tens of millions 
| Norris et al.|[20lT |, within such surveys we expect to find 
between tens to hundreds of thousands of Head-Tailed galax- 
ies \Mao et al.|2011|[Dehghan et al.|2011[ ), which will form 
an important population for cluster physics. In recent years, 
the use of HTs as probes of local over densities in their near 
by environment has increased somewhat. This includes as 
signposts for detecting clusters ( Blanton et al.|2003 Smolcic 
|2007[ ) and local over densities ( Mao et al.||2010| ). However, 
despite this oft cited use, few studies have fully utilized HTs 
as environmental probes. This work highlights how rich HT 
galaxies will be for understanding clusters in next generation 
radio surveys. 
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